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ABSTRACT 

We investigate whether it is possible to explain the wide range of observed gradients 
in early type galaxies in the framework of monolithic models. To do so, we extend 
the set of hydrodynamical simulations by Pipino et al. (2008a) by including low- 
mass ellipticals and spiral (true) bulges. These models satisfy the mass-metallicity and 
the mass-[Q;/Fe] relations. The typical metallicity gradients predicted by our models 
have a slope of -0.3 dex per decade variation in radius, consistent with the mean 
values of several observational samples. However, we also find a few quite massive 
galaxies in which this slope is -0.5 dex per decade, in agreement with some recent data. 
In particular, we find a mild dependence from the mass tracers when we transform 
the stellar abundance gradients into radial variations of the Mg2 line-strength index, 
but not in the Mgi^. We conclude that, rather than a mass- slope relation, is more 
appropriate to speak of an increase in the scatter of the gradient slope with the galactic 
mass. We can explain such a behaviour with different efficiencies of star formation in 
the framework of the revised monolithic formation scenario, hence the scatter in the 
observed gradients should not be used as an evidence of the need of mergers. Indeed, 
model galaxies that exhibit the steepest gradient slopes are preferentially those with 
the highest star formation efficiency at that given mass. 



1 INTRODUCTION 



Negative metallicity radial gradients in the stellar popula- 
tions are a common feature in spheroids (e.g. Carollo et al. 
1993, Davies et al. 1993) and must be predicted by every 
theory for the formation of elliptical galaxies. A possible fin- 
gerprint of a given galaxy formation scenario might be the 
(lack of) correlation between gradient properties (e.g. the 
slope) and either global galactic properties (namely mass, 
stellar velocity dispersion a, total magnitude) or central ones 
(e.g. central metallicity or [< q/Fo]). From the theoreti- 
cal point of view, in fact, steep metallicity gradients are ex- 
pected from classical dissipative collapse models (e.g. Larson 
1974, Chiosi & Carraro, 2002) and their [revised) up-to-date 
versions which start from semi- cosmological initial condi- 
tions (e.g. Kawata 2001, Kobayashi, 2004). The abundance 
gradient arises because the stars form everywhere in a col- 
lapsing cloud and then remain in orbit with a little inward 
motion|j whereas the gas sinks further in because of dis- 
sipation. This sinking gas contains the new metals ejected 
by evolving stars so that an abundance gradient develops in 
the gas. As stars continue to form their composition reflect 

^ Stars will spend most of their time near the apoccntrc of their 
orbit. 



the gaseous abundance gradient. The original dissipative 
models predict a steepening of the gradient as the galac- 
tic mass increases, mainly because the central metallicity is 
quickly increasing with mass Q whereas the global one has 
a milder variation (Carlberg 1984). At the same time they 
predict metallicity gradient as steep as -0.5 dex per decade 
variation in radius. On the other hand, the few attempts 
to study the gradients in the merger-based models hint for 
very shallow (if any) gradient (Bekki & Shyoia, 1999), less 
steep than the mean observational values and than the pre- 
dictions from monolithic collapse models. Moreover, it seems 
that dry mergers flatten pre-existing gradients (Di Matteo 
et al., 2009). Indeed, when the two scenarios (monolithic 
collapse and mergers) are considered as two possible chan- 
nels working at the same time, the scatter in the predicted 
gradients for such a population of galaxies seem to be in 
agreement with observations (Kobayashi, 2004). 

More recently, observations showed that successful 
models for elliptical galaxies should also reproduce the [< 
a/Fe>]-mass relation (Worthey et al. 1992, Thomas et al. 



■^ The fit of the mass-metallicity relation, namely the increase of 
the mean metal content in the stars as a function of galactic mass 
(O'Connel, 1976), was the main success of these original models. 
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2007) as well as the observed gradients in the [< a/Fe >] 
ratios (Mehlert et al. 2003, Annibali et al. 2007, Sanchez- 
Blazquez et al. 2007, Rawle et al., 2008). Indeed, these ob- 
servations show that the slope in the [< a/Fe >] gradient 
has a typical value close to zero and does not correlate with 
mass. 

These observations have been interpreted by Pipino et 
al. (2008a; Paper I, hereafter) ID hydrodynamical code as 
the fact that the suggested outside-in mechanism for the 
formation of the ellipticals is not the only process respon- 
sible for the formation of gradients in the abundance ra- 
tios. Other processes should be considered such as the in- 
terplay between the star formation (SF) timescale and gas 
flows. While such an interplay flattens the [< a/Fe >] gra- 
dient to the value required by observations, it still enables 
galaxies to harbor gradients in [<Fe/H>] and [<Z/H>] in 
agreement with the most recent observations (see Section 2). 
Pipino et al. (2008b; Paper II) calibrated such a model by 
means of the resolved stellar populations in the Milky Way 
bulge. As a matter of fact, spiral iruqj bulges remarkably 
follow many fundamental constraints for ellipticals such as 
the mass-metallicity and the mass-[< a/Fe>] relations (see 
below), the only difference being that they might be rejuve- 
nated systems (Thomas & Davies, 2006). 

The aim of this paper is to explore a wider range of cases 
by extending the analysis of Paper I to lower masses, includ- 
ing bulges, and compare them to the latest observational 
results. In this way we can study the correlation between 
gradient slopes and galactic mass (if any) in order to under- 
stand whether the monolithic galaxy formation scenario is 
in agreement with the recent observational evidences. 

In Sec. [2] we give a brief overview of the observations 
regarding metallicity gradients in ellipticals. The main char- 
acteristic of the model are briefly described in Section [S] We 
characterize the global properties of our models in Section 
4, present our results in Section [5] discuss them in Section 
6 and draw our conclusions in Section [T] 



2 THE OBSERVATIONAL BACKGROUND 

In general, observations show that the majority of ellip- 
ticals has as typical decrease in metallicity of 0.2-0.3 dex 
per decade in radius (e.g. CaroUo et al. 1993, Davies et al. 
1993). However, a large scatter in the gradient slope at a 
given galactic mass is also observed. The exact slope de- 
pends on the line-strength index used to infer the metal- 
licity. Below we give a brief historical perspective for what 
concerns the relation between gradient slope and mass. We 
refer the reader to other works (e.g. Sanchez-Blazquez et al. 
2006) for a review about the debate on the observations in 
the literature. 

Indeed, a positive correlation of the metallicity gradi- 
ent slope with the galactic mass - namely gradients becom- 
ing more negative at higher galactic masses - (in agreement 
with Larson 1974's prediction), has been reported by Car- 
oUo et al. (1993), but only for masses lower than 10^ ^Mq. 
In fact, CaroUo et al. (1993) found a Hattening of the ob- 
served gradients in the most massive galaxies of their sample 

^ In the rest of the paper we will consider only the class of true 
bulges, (Korniendy & Kennicutt, 2004). 



and ascribed this fact to: i) an increase in the importance 
of mergers; or ii) a less important role of dissipation in the 
formation of the most massive galaxies. The positive correla- 
tion of the slope with the galactic mass was later confirmed 
by some authors (e.g. Gonzalez & Gorgas, 1996) over the 
entire mass range and denied by others who either found no 
statistical evidence for such a correlation (e.g. Kobayashi & 
Arimoto, 1999) or a very mild opposite trend (e.g. Annibali 
et al., 2007). We notice that several of the studied sam- 
ples were quite small or not homogeneous (e.g. Kobayashi 
& Arimoto, 1999). In recent years, a positive correlation of 
gradient slope with mass has been suggested again by Forbes 
et al. (2005), Sanchez-Blazquez et al (2007), for the entire 
mass range of elliptical galaxies. Ogando et al. (2005), rather 
than a clear trend, noticed an increasing number of E and 
SO galaxies harboring steep Mg2 gradients with increasing 
velocity dispersion. Interestingly, Spolaor et al. (2009) found 
a similar result for massive ellipticals, whereas, for the first 
time, detected a clear gradient slope-mass relation at low 
mass end (Fornax and Virgo dwarf). Spolaor et al's re- 
sult has been questioned by Koleva et al. (2009a), who do 
not observe any such a trend in another sample of dwarf 
galaxies in the Fornax cluster. To date, no one has offered a 
convincing explanation for the discrepancy between obser- 
vational results (unfortunately Koleva et al.'s and Spolaor 
et al's samples do not overlap!). One problem, of course, 
is the small number statistics. Issues related to the reduc- 
tion and analysis process have been excluded as a cause for 
this discrepancy (Koleva et al., 2009b). Moreover, as we will 
discuss later in the comparison between our models and ob- 
servations, different authors use different (combinations of) 
indices to estimate the age, a/Fe and metallicity indices. 
This is sufficient to maJie the inferred gradients appear either 
stronger or weaker (e.g. Sanchez-Blazquez et al., 2006). In 
addition, they use different SSP libraries and minimization 
techniques to transform their data into metallicity (either 
[Z/H] or [Fe/H]) and ages, thus introducing further issues in 
the interpretation (see Pipino et al., 2006 for an extended 
discussion). 

Bulges have gradients in metallicity (Goudfroij et al. 
1999, Proctor et al. 2000) and [< a/Fe >] ratios (Jablonka 
et al. 2007) with the same properties as those in ellipticals. In 
particular, Jablonka et al. (2007) described the variation in 
the gradient slope as a function of mass as a multi-step pro- 
cess rather than a smooth transition in gradient amplitude 
with velocity dispersion. According to the latter authors, 
at large masses the dispersion among gradients is large but 
small gradients are relatively rare. At smaller masses, in- 
stead, galaxies with very weak gradients appear in larger 
number. 



3 THE MODEL 

We adopted a one-dimensional hydrodynamical model 
(Frankenstein) that follows the time evolution of the den- 
sity of mass (p), momentum (m) and internal energy (e) 
of a galaxy, under the assumption of spherical symmetry. 
In order to solve the equation of hydrodynamics with a 
source term we made use of the code presented in Ciotti 
et al. (1991), which is an improved version of the Bedogni 
& D'Ercole (1986) Eulerian, second-order, upwind Integra- 
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tion scheme (see their Appendix). Here we report the gas- 
dynamics equations: 



^ + V • [pu) = ap, - *, 



+ V ■ {q^u) = a'p* - ^d' Ip, 



dt 

dm „ , , , , , „ 

— + V • [mu) = P5 - (7 - 1)V£ - *«, 



(1) 
(2) 
(3) 



-^ + V-{£u) = -(7-l)eV-M-L+ap, i ^o + -u j -*e/p.(4) 

The parameter 7 = 5/3 is the ratio of the specific heats, g 
and It are the gravitational acceleration due to the total 
mass distribution (stars and dark halos) and the fluid veloc- 
ity, respectively. The source terms on the r.h.s. of equations 
(l)-(4) describe the injection of total mass and energy in the 
gas due to the mass return and energy input from the stars. 
a{t) — Q,(i) + QsNii(i) + asNia(i) is the sum of the spe- 
cific mass return rates from low-mass stars and SNe of both 
Type II and la, respectively, eo = SkTo / {2pmp) is the in- 
jection energy per unit mass due to the SN explosions, and 
To is the injection temperature. The positive source term 
on the right-hand side of the energy equation describes the 
heating of the gas by SN blast waves and by the relative 
motion of the mass-losing stars and the ISM (kinetic heat- 
ing). ^ is the astration term due to SF. L — nenpA(T, Z) is 
the cooling rate per unit volume, where for the cooling law, 
A{T,Z), we adopt the Sutherland & Dopita (1993) curves. 
This treatment allows us to implement a self-consistent de- 
pendence of the cooling curve on the metallicity (Z) in the 
present code. We do not allow the gas temperature to drop 
below 10* K, as the Sutherland & Dopita (1993) functions 
are calculated only above this limit. We are aware that fix- 
ing the minimum gas temperature can be a limitation of the 
model, but this is done in order to avoid the complexity of 
the cooling at lower temperatures. Moreover, as it can be 
seen from Paper I Figs. 1 and 2, at the time of the occur- 
rence of the winds (and actually for most of the pre- wind 
evolution) the majority of the models exhibit T >> 10'* K. 
£)' represents the mass density of the i — th element, 
and a* the specific mass return rate for the same element. 



where tcooi and t// are the local cooling and free-fall 
timescales, respectively, and esp is a suitable SF parame- 
ter that contains all the uncertainties on the timescales of 
the SF process that cannot be taken into account in the 
present modelling and will be taken as a free parameter in 
our models. In fact, star formation is an inherently 3D pro- 
cess which cannot be even approximately simulated by ID 
simulations. Moreover, star formation occurs on small scale, 
much smaller than any possible mesh resolution when the 
whole galaxy must be covered by the numerical grid. We re- 
call that the final efficiency, namely the fraction of gas that 
eventually turned into stars, is an output of the model. 

We assume that the stars do not move from the grid- 
points at which they have been formed, since we expect that 
the stars will spend most of their time close to their apoc- 
entre. 

3.1 Chemical evolution 

The nucleosynthetic products enter the mass conservation 
equations via several source terms, according to their stel- 
lar origin. A Salpeter (1955) initial mass function (IMF) 
constant in time in the range 0.1 — 50Mq is assumed. We 
adopted the yields from Iwamoto et al. (1999, and references 
therein) for both SNIa and SNII. The SNIa rate for a SSP 
formed at a given radius is calculated assuming the single 
degenerate scenario and the Matteucci & Recchi (2001) de- 
lay time distribution (DTD). These quantities, as well as 
the evolution of single low and intermediate mass stars, were 
evaluated by adopting the stellar lifetimes given by Padovani 
& Matteucci (1993). The solar abundances - used to present 
our values in the "[< >]"notation - are taken from Asplund 
et al. (2005), unless otherwise stated. Note that, as far as 
gradient slope are concerned, the actual solar scale does not 
make any difference. 

In order to study the mean properties of the stellar com- 
ponent in ellipticals, we need average quantities related to 
the mean abundance pattern of the stars, which, in turn, can 
allow a comparison with the observed integrated spectra. 
In particular, we make use of the luminosity-weighted mean 
stellar abundances. Following Arimoto & Yoshii (1987), we 
have; 



with J^Zi "' = «• Eq- (2) represents a subsystem of four < O/^^' >v= 2^ nkj(0/Fe)iLv,k/ 2^ Uk.iL^ 



(6) 



equations that follow the hydrodynamical evolution of four 
different ejected elements (namely H, He, O and Fe). This 
set of elements is good enough to characterize our simulated 
elliptical galaxy from the chemical evolution point of view. 
We divide the grid in 550 zones 10 pc wide in the inner- 
most regions, and then slightly increasing with a size ratio 
between adjacent zones equal to 1.03. At the same time, 
however, the size of the simulated box is roughly a factor of 
10 larger than the stellar tidal radius. This is necessary to 
avoid possible perturbations at the boundary affecting the 
galaxy and because we want to have a surrounding medium 
that acts as a gas reservoir for the models. We adopted a 
reflecting boundary condition in the center of the grid and 
allowed for an outflow condition in the outermost point. 

At every point of the mesh we allow the SF to occur at 
the following rate: 



*: 



up- 



<^SF 



max {t cool, tff) 



(5) 



where Uk^i is the number of stars binned in the interval 
centered around (0/Fe)i with V-band luminosity Lv,k. We 
then take the logarithm and express the quantities in so- 
lar units. Similar equations hold for [< Fe/H >v] and the 
global metallicity [< Z/H >v]- Generally the mass averaged 
[Fe/H] and [Z/H] are slightly larger than the luminosity av- 
eraged ones, except for large galaxies (see Yoshii & Arimoto, 
1987). We will present our results in terms of [< Fe/H >v] 
and [< Z/H >v], because the luminosity- weighted mean 
is much closer to the actual observations and might differ 
from the mass-average, unless otherwise stated. Therefore 
we drop the subscript V in the remainder of the paper. 

3.2 Model classification and initial conditions 

The initial set-up of the new simulations for low-mass ellip- 
ticals is presented in Table 1, where the name of the model. 
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the gas density {pcore,gas ) as well as the initial gas tempera- 
ture, the star formation parameter esp and the Dark Matter 
(DM) halo mass are reported. In the same table we include 
also the models already presented in Paper I and the bulges 
(see below). 

We recall that in Paper I we defined the following two 
families of models according to the total initial DM and gas 
content: Model M - 2.2 • lO^^M© DM halo and ~ 2 ■ lO^Af© 
of gas; Model L - 5.7 ■ 10^^ M© DM halo and ~ 6.4 • lO^M© 
of gas. The DM potential has been evaluated by assuming 
a distribution inversely proportional to the square of the ra- 
dius at large distances (Silich & Tenorio-Tagle 1998). These 
quantities have been chosen to ensure a final ratio between 
the mass of baryons in stars and the mass of the DM halo 
of around 0.1. 

In order to model ellipticals less luminous than the ones 
presented above, we assume that the galaxy assembly oc- 
curs in a 0.3 times smaller and 0.1 times lighter DM halo. 
This guarantees that we model ~ 0.5 — 2 x 10^° M© galaxies 
(stellar mass) in ~ 2 x 10^^ M© Dark Matter haloes. Note, 
however, that the final mass in stars is determined by the 
interplay between the SF efficiency and the duration of the 
SF process itself (regulated by infall and stellar feedback). 
Therefore we may have two galaxies with the same DM po- 
tential and differing stellar masses because of the different 
evolutionary paths. 

Concerning the bulges, instead, we assume that they 
are stellar systems with mass ~ 2 x lO^'^M© embedded in 
a ~ 100 times more massive Dark Matter halo, since bulges 
occupy only the central part of their large hosts. We ne- 
glect the presence of a disc, which requires a much longer 
timescale to be built (e.g. Zoccali et al., 2006 from the ob- 
servational viewpoint, Matteucci & Brocato, 1990, Ballero 
et al.,20Q7, from the theoretical one). Moreover, Sarajedini 
& Jablonka (2005) suggest a common scenario for the for- 
mation of bulges that is not linked to the host galaxy forma- 
tion. Finally, the observations we are comparing our results 
to have been derived by accurately selecting edge-on galax- 
ies. Therefore the contamination from the discs should be 
minimal. 

To generate different models we mainly vary the gas 
temperature and the efficiency of star formation as well as 
the initial gas density distribution. 

In particular, the gas can initially be an isothermal 
sphere (models fiagged as IS) in equilibrium within the 
galactic potential well (i.e. due to both DM and gas). The 
actual initial temperature is lower than the virial tempera- 
ture, in order to induce the gas to collapse. This is an ex- 
treme case in which we let all the gas be accreted before 
the SF starts. In other models, instead, the gas has uniform 
distribution within the whole computational box (models 
fiagged as flat). At variance with the previous models, in 
this case we let the SF process start at the same time as 



the gas accretion. The values for pc 



are set in order 



not to have too much gas in the grid, namely higher than 
the typical baryon fraction in high density environment (i.e. 
1/5-1/10 as in galaxy cluster, e.g. McCarthy et al. 07). 

The initial gas temperature ranges from 10**^^ K (cold- 
warm gas) to 10®^^ K (virialised haloes). This range of 
temperature is consistent with the typical findings of simu- 
lations of high-redshift galaxy formation. In fact, a common 
assumption in galaxy formation models has that the gas ac- 



creted by a DM halo is shock-heated to the host halo virial 
temperature (lO^K), and only then is able to cool down and 
feed star formation (e.g. White & Rees 1978). This scenario 
justifies the models with a high initial temperature and IS 
gas profile, with the only difference that the amount gas 
reservoir is not regulated by any "cosmological" infall his- 
tory. Slightly different (high) initial temperatures may be 
used to regulate/delay the infall rate on the actual pro- 
togalaxy. We note that the gas cools very rapidly, there- 
fore the actual starting value is less important than, e.g., 
the chosen esp- Recent simulations show that the gas may 
be accreted through cold filaments (e.g. Dekel & Birnboim, 
2006) streaming through the shock-heated gas. In this case 
it will have temperatures of about 10'*~^K. The majority of 
the models presented here {flat profile and warm tempera- 
ture) are motivated by these recent results. The ID nature of 
our study hampers us to mode these "cold" accretion fiows, 
therefore we simply varied intial gas density and temper- 
ature in order to give a reasonable approximation to this 
picture. 

In general, we assume values for the star formation 
parameter between 0.1 and 10. These values guarantee star 
formation rates of 10-500 M©/yr (c.f Paper I, Fig. 8) in mas- 
sive galaxies, comparable with the observations of high red- 
shift star forming objects. A preliminary exploration of the 
parameter space returned that smaller values of the star for- 
mation parameter give rise to too extended star formation 
histories (and hence too low [a/Fe] ratios). On the other 
hand, higher values would lead to too small (in terms of 
stellar to total mass ratio) and too large (in terms of Reff) 
galaxies - similarly to a what happens when one adopts a 
100% SN efficiency (cf model MaSN, Paper I). In fact, the 
strong feedback from SNe halts the SF too early by prevent- 
ing further accretion of gas. Such a galaxy would also have 
a too high [a/Fe] . These models have been discarded during 
the preliminary analysis that led to Paper I. Both the SNIa 
and SNII efficiency is assumed to be constant esN = 0.1 (see 
Paper I, Pipino et al. 2005). 

We choose Reff,* as the radius that contains 1/2 of the 
stellar mass and, therefore, is directly comparable with the 
observed effective radius, whereas we will refer to Rcore,* as 
the radius encompassing 1/10 of the galactic stellar mass. 
We did not fix Rcore,* = O.lReff,* a priori, in order to have 
a more meaningful quantity, which may carry information 
on the actual simulated stellar profile. In most cases, this 
radius will correspond to '^ 0.05 — 0.2Rsff,*, which is the 
typical size of the aperture used in many observational works 
to measure the abundances in the innermost regions of el- 
lipticals. 

Finally, we did use the following notation for the metal- 
licity gradients in stars Aq/pb ~ ([< O/Fe >]core — [< 
O/Fe >]eff)/log{Rcore,»/Reff,t)', a similar expression ap- 
plies for both the [< Fe/H >] and the [< Z/H >] ratios. 
Hence, the slope is calculated by a linear regression between 
the core and the half-mass radius, unless otherwise stated. 
Clearly, deviations from linearity can affect the actual slope 
at intermediate radii (see Fig. 3). 

For all the models the velocity dispersion a is evaluated 
from the relation M = 4.65 ■ 10^ (cr/kms"^)^ Reff/kpcMe 
(Burstein et al., 1997). We warn the reader that we assume 
that our model galaxies are virialised objects in order to 
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Table 1. Input parameters 



(10- 



gcm ) 



initial 
profile 



(K) 



Massive ellipticals (Paper I) 



Mai 

Ma2 

Ma3 

Mbl 

Mb2 

Mb3 

Mb4 

La 

Lb 



0.6 


IS 


1 


10'' 


22 


0.6 


IS 


10 


10" 


22 


0.6 


IS 


2 


10" 


22 


0.06 


flat 


1 


10' 


22 


0.2 


flat 


1 


10^ 


22 


0.06 


flat 


10 


10*^ 


22 


0.6 


flat 


1 


lO" 


22 


0.6 


IS 


10 


10^ 


57 


0.6 


flat 


10 


lO" 


57 


0.3 


flat 


0.5 


10^ 


2 


0.3 


IS 


0.5 


10^ 


2 


0.3 


IS 


3 


10^ 


2 


0.01 


flat 


ID'S 


2 




0.03 


flat 


1 


10^ 


2 


0.02 


flat 


1 


10^ 


2 


0.02 


flat 


0.1 


10^ 


2 


0.02 


flat 


0.3 


10^ 


2 


0.02 


flat 


0.2 


10^ 


2 


0.02 


flat 


0.3 


lO** 


2 


0.02 


flat 


0.2 


lO" 


2 


0.02 


flat 


1 


lO" 


2 


0.007 


flat 


1 


10^ 


2 


0.007 


flat 


0.1 


lO'"- 


2 


0.02 


flat 


3 


10^ 


2 


0.02 


flat 


10 


10^ 


2 



Low mass ellipticals 



Ela 

Elb 

Elc 

E2a 

E2b 

E2c 

E2d 

E3a 

E3b 

E3c 

E3d 

E4a 

E4b 

E4c 

E5 

E6 



Bulges 



bulgel 
bulge2 
bulge3 
bulge4 
bulge5 



0.02 


IS 


1 


10" 


20 


6. 


IS 


3 


10^ 


20 


6. 


IS 


3 


lO" 


20 


0.02 


flat 


3 


105 


20 


0.007 


flat 


3 


10^ 


20 



Models called E are low-mass ellipticals, whereas models called bulge are spiral bulges. The flags flat and 
IS pertain to the initial gas distribution ^vhich can be either constant \vith radius (flat) or an isothermal 
sphere (IS), respectively. The model hulgeS has been used in Paper II for a calibration on the chemical 
properties of the resolved stars in the Milky Way bulge. 



assign them a stellar velocity dispersion from their mass and 
effective radius, because we do not model stellar kinematics. 



RESULTS: 
MODELS 



GLOBAL PROPERTIES OF THE 



We start the analysis of our results by briefly discussing 
some general properties that hold for the entire sample of 
models - i.e. ellipticals and bulges - whose relevant predicted 
properties are listed in Table [2] (including massive ellipticals 
from Paper I). In particular, we show the final (i.e. after SF 
stops) values for the stellar mass and effective radius, the 
[< O jFe >] abundance ratio in the galactic centre and the 
gradients in [< OjFe >] and [< ZjH >]. 

The relation between [a/Fe] and mass tracers (see e.g. 
Worthey et al. 1992, Nelan et al. 2005, Thomas & Davies, 
2006), is satisfied, as shown in Fig. [T] It is important to 
ensure that the models fulfill such a relation, as it is the 
most severe test-bench for a galaxy formation scenario (see 
Pipino et al. 2009a). We note that the mass-metallicity re- 
lation is also satisfied, since our massive objects have an 
average stellar metallicity which is super-solar, whereas the 
simulated low mass ellipticals and bulges have solar metal- 



licity at most. More quantitatively, a linear fit to our model 
predictions gives [<Z/H>]core = — 1.14 -f 0.57/o(;cr to be 
compared with the relation [<Z/H>]core = — 1.06+0. SSZop cr 
inferred by Thomas et al. (2005) within the same aperture 
for observed ellipticals. The robustness of our predictions is 
supported by the fact that our models obey to the above 
mentioned observational constraints. This ensures that we 
investigate the relation between abundance and abundance 
ratios gradients by means of models that are able to re- 
produce the main chemical properties of the ellipticals. Re- 
markably, the above-mentioned relations are in place already 
after 0.5 - 1 Gyr since the beginning of the star formation. 
In the two following sections we highlight other main 
features of model ellipticals and bulges, respectively. 



4.1 Elliptical galaxies 

In brief, we first recall from Paper I how the formation of 
a galaxy proceeds in our model. We take the case La as an 
example. At times earlier than 300 Myr the gas is still accu- 
mulating in the central regions where the density increases 
by several orders of magnitude, with a uniform speed across 
the galaxy. The temperature drops due to cooling, and the 
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Figure 1. The [< O/ Fe >\core-o' relation predicted for our model 
ellipticals and bulges shown as crosses in the figure. Data from 
Thomas et al. (2007) are shown as contours. A liner regression to 
Thomas et al. (2007) is shown by a thin solid line. Note that in this 
plot we rescaled our [< O/ Fe >] values in order to be consistent 
with the solar abundances used by Thomas et al. (2007). Spiral 
bulges obey to the same relation (Thomas & Davics, 2006). 



Figure 2. Metallicity versus v^sc gradient slope for our models. 
The central Mgi, value for each model galaxy is given by an aster- 
isk, whereas the value at 1 i?e// by a cross. Each couple of points 
connected by a line represents a galaxy (the local relation). The 
dashed line is the global relation (Scott et al., 2009) along with 
the Scr dispersion (dotted lines). 



SF can proceed at a very high rate (~ lO^~^M0yr~^), at 
variance with the outermost regions, that complete their 
build-up in the first 100 Myr. This implies that a metal 
rich medium, dominated by SNIa ejecta, pollutes the gas 
supply for the SF in the inner regions. After 400 Myr, the 
gas speed becomes positive (i.e. out-flowing gas) at large 
radii, and at 500 Myr almost the entire galaxy experiences 
a galactic wind. At roughly 1.2 Gyr, the amount of gas left 
inside the galaxy is below 2% of the stellar mass. This gas 
is very hot (around 1 keV) and still flowing outside. 

The galactic wind occurs first in the outer regions and 
then in the more inner zones of the galaxy because the work 
to extract the gas from the outskirts is less than the work 
to extract the gas from the center of the galaxy. The age 
differences between internal and external zones, however, 
are less than 1 Gyr and this ensures that our models are 
globally Q-enhanced. In this way our models are consistent 
with the observed age gradients (references in Sec. 2) and 
with the [a/Fe]-mass relation. The picture sketched above 
applies to the lower mass models presented here. The fact 
that in our galaxy formation scenario the metallicity gra- 
dients arise because of the different times of occurrence of 
galactic winds in different galactic regions implies that the 
stellar metallicity is a function of the local escape velocity 
Vesc for all the galaxies. In fact, in the regions where Vesc 
is low (i.e. where the local potential is weaker), the galactic 
wind develops earlier and the gas is less processed than in 
the regions where v^sc is higher (see Martinelli et al. 1998). 
Such a relation as been originally suggested by several au- 
thors (e.g. Davies et al., 1993, Peletier et al., 1990) and now 
confirmed by Scott et al. (2009). Here we can also show that 
the local index-Uesc trend matches the global scaling (Scott 
et al., 2009). In particular, we make use of the definition 
Vssc ~ •\/(— 2$), where $ is the potential due to stars and 
DM, in agreement with the definition used by the observers. 
Some important caveats apply to this comparison. In ob- 
servations, Vesc depends on the modelling of the potential. 



Moreover, our models are spherically symmetric, whereas 
observed galaxies are not. In Fig. [2] we show that metallic- 
ity (given by the index Mgb) versus Vesc gradient slope for 
our models. The central Mgt, value for each model galaxy 
is given by an asterisk, whereas the value at 1 Reff by a 
cross. Each couple of points connected by a line represents 
a galaxy: this is the local relation. The dashed line is the 
observational global (i.e. the fit to the central values of Mgt, 
and Vesc in observed galaxies) trend reported by Scott et al. 
(2009) along with the "ia dispersion (dotted lines) . We show 
that the models presented in this paper reproduce the ob- 
served trend within the observed scatter. The fact that the 
each galactic region follows the global trend strongly sug- 
gests the idea that a uniform process - like the monolithic 
collapse - is behind the formation of the gradients. 



4.2 Galaxy bulges 

Remarkably, all the results discussed in the previous sections 
apply to smaller objects (but embedded in much more mas- 
sive haloes) such as the galaxy bulges, although the gradient 
slopes are slightly smaller (see entries in Table [Sjl. The main 
difference is that, due to their host galaxy potential well, 
strong and long lasting winds do not develop. We also find 
that the bulge formation is fast in agreement with the orig- 
inal suggestion by Matteucci & Brocato (1990), Elmegreen 
(1999) and the more recent work by Ballero et al. (2007). 

We take advantage of the classical bulges as a further 
tool to calibrate our models. Indeed, in Paper II (where we 
refer the reader for further details) we compared our model 
predictions to the properties of the resolved stellar popula- 
tion observed in the Milky Way bulge by using the model 
bulges and found a remarkable agreement. This model has a 
stellar mass of ~ 2 x 10^" Mq and a radius of ~ Ikpc in order 
to match the observed properties of our own Galaxy bulge 
(e.g. Minniti & Zoccali, 2008). The same model reproduces 
the chemical constraints coming from the Bulge integrated 
light, in that it predicts the following values for the indices 
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Table 2. Model results 



Model 


(lOlOAfg) 


(kpc) 


[< O/Fe >.,cor, 


= 1 ^O/F^ 


^Z/H 


Massive ellipticals (Paper I) 


Mai 


6.0 


12 


0.29 


0.02 


-0.19 


Ma2 


25. 


7.7 


0.22 


-0.21 


-0.52 


Ma3 


26. 


8.3 


0.3S 


-0.17 


-0.03 


Mbl 


6.0 


17 


0.14 


0.09 


-0.20 


Mb2 


3.0 


8.7 


0.33 


0. 


-0.18 


Mb3 


21 


8.8 


0.17 


-0.08 


-0.34 


Mb4 


26 


5.4 


0.42 


-0.08 


-0.20 


La 


26 


29 


0.14 


0.19 


-0.50 


Lb 


29 


21 


0.12 


0.32 


-0.30 


Low mass ellipticals 


Ela 


0.74 


1.7 


0.08 


-0.04 


-0.26 


Elb 


0.74 


1.7 


0.36 


-0.13 


-0.21 


Elc 


0.74 


1.7 


0.28 


-0.11 


-0.21 


E2a 


1.5 


0.9 


0.19 


-0.03 


-0.29 


E2b 


1.8 


0.6 


0.14 


-FO.Ol 


-0.29 


E2c 


1.4 


0.89 


0.26 


-0.04 


-0.30 


E2d 


0.27 


2.3 


0.17 


-1-0.01 


-0.29 


E3a 


0.88 


1.6 


0.18 


-0.005 


-0.27 


E3b 


0.65 


1.1 


0.11 


-1-0.07 


-0.28 


E3c 


0.93 


1.6 


0.09 


-0.01 


-0.32 


E3d 


0.6 


1.1 


0.03 


-HO. 06 


-0.25 


E4a 


1 


1.7 


0.16 


-0.21 


-0.33 


E4b 


0.35 


0.6 


0.22 


-0.04 


-0.36 


E4c 


0.05 


0.5 


0.17 


-0.01 


-0.22 


E5 


1 


1.7 


0.16 


-0.20 


-0.38 


E6 


1 


1.7 


0.11 


-0.10 


-0.34 


Bulges 


bulge 1 


0.06 


2 


0. 


0.09 


-0.36 


bulge2 


1.8 


1 


0.40 


-0.07 


-0.22 


bulge3 


2.3 


0.8 


0.3 


0.07 


-0.36 


bulge4 


3.7 


0.7 


0.29 


0.00 


-0.37 


bulges 


1.0 


0.4 


0.2S 


0.00 


-0.30 



Models called E arc low-mass ellipticals, whereas models called bulge are spiral bulges. Values predicted 
after the SF has finished. 



H^ = 1.61, Mg2=0.29 and < Fe >=2.46 in good agreement 
with the observed values of H^ = 1.5±0.6, Mg2=0.23±0.04 
and < Fe >=2.15±0.4 (Puzia et al., 2002). This is an im- 
portant point that must be stressed: the abundance (and 
abundance gradients) that may be inferred from the analysis 
of Lick indices are average values. With resolved stellar pop- 
ulations (Paper II) is possible to show that the models pre- 
sented here not only explains the average values, namely the 
mean properties of a Composite Stellar Population (CSP), 
but also their evolution in the [0/Fe]-[Fe/H] plane, namely 
the composition of each SSPs that make a CSP. Our fidu- 
cial model assumes Salpeter (1955) IMF, which successfully 
reproduces the properties of massive spheroids. In Paper II 
we show that the stellar metallicity distribution predicted 
by such a model reproduces the observed K-giant metallicity 
distributions for the Milky- Way bulge. We refer to Paper II 
(c.f. Fig. 2) for the test of other possible IMFs, motivated by 
either observations or theoretical efforts, which seems more 
appropriate for bulges. 



The reader should note that we present several other 
models for bulges which do not necessarily have properties - 
such as stellar mass or radius - similar to those of the Milky 
Way bulge. 



5 THE PREDICTED GRADIENTS 

In this section we turn our attention on the main topic of 
this work and investigate the possible dependence of the gra- 
dient slope - and its scatter - from either the stellar mass 
or some mass tracers. We first focus on the actual predic- 
tion from the modeller's point of view, namely gradients 
in abundance and abundance ratios versus mass and central 
[<0/Fe>], whereas we refer to Sec. 5.2 for our model predic- 
tions transformed into observational line-strength indices. 

The metallicity profiles predicted by our model ellipti- 
cals over the 0.1-lReff range are shown in Fig. [S] In the 
vast majority of the cases we predict a linear decrease of 
the metallicity with log (r) and thus justifies the adopted 
definition for Az/h- We refer the reader to Paper I (Fig. 
8) comparison between the observed and predicted [O/Fe] 
profiles for some relevant cases. 

For elliptical galaxies we make use of Mehlert et al.'s 
(2003) and Annibali et al.'s (2007) datasets, whose samples 
are larger than Ogando et al.'s one, although the former do 
not find such a strong correlation between gradient slope 
and mass as the latter (other works with less galaxies are 
not taken into account in order not to have a poor statis- 
tics). For bulges we adopt the data from Jablonka et al. 
(2007), who explore a range in velocity dispersions similar 
to the above- mentioned articles. Unfortunately, we cannot 
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Figure 3. Mctallicity profiles predicted by our model ellipticals. 



use a homogeneous set of observables to constrain both the 
theoretical and the observational predictions for several rea- 
sons. In the first place, in several articles the authors do not 
extract the [<0/Fe>] abundance ratio gradient from their 
line-strength indices (e.g. Ogando et al. 200qJ, Kobayashi & 
Arimoto 1999, and Sanchez-Blazquez et al. 2006). Secondly, 
in all cases the stellar mass is not observed, whereas only the 
stellar velocity dispersion is given. Finally, several authors 
rely on a different sub-set of the Lick line-strength indices 
to infer the metallicity. 



5.1 Theoretical relations with mass and mass 
tracers 

With the above mentioned caveats in mind, in Fig. |4] we 
present our predictions regarding the theoretical relation be- 
tween abundance gradients and mass tracers (namely the 
stellar velocity dispersion and the central [< OjFe >]). The 
remainder of this section is devoted to fully describe Fig. |4l 



5.1.1 Gradients in metallicity 

Let us first focus on the upper row of Fig. [l] the total metal- 
licity gradient. In the left panel we show the [< Z/H >] 
gradient slope in the stellar component predicted by our 
model for ellipticals (hollow circles) and bulges (full dots) 
as a function of the stellar mass when all galaxies are con- 
sidered. This is the actual prediction of our models. Formal 



* Noticeably, they could not convert the indices into abundances 
in several galaxies whose combination of index values fell outside 
Thomas et al. (2003) SSP libraries. We refer the reader to Pipino 
et al. (2006) and Paper I for a detailed discussion on the theo- 
retical aspects of such a problem. Here we just mention that SSP 
libraries do not cover all the possible combinations in the space 
[< O/Fe >]- [< Fe/H >]-[< Z/H >], being typically built just 
as functions of two of them. 



linear regression fits to the entire sample of model galaxies 
(solid line), to the galaxies with steepest gradients (dotted 
line ) and to dwarf ellipticals (dashed line) are shown. We 
predict a very mild trend in mass. In the high mass region, 
our model predictions span a range in the gradient slopes 
similar to the observed values. Neither our models, nor the 
three observational samples (taken togheter) show any sign 
of (anti-)correlation as suggested by the single authors. We 
therefore conclude that, in this mass range, is more appro- 
priate to speak of an increase in the scatter of the gradient 
slope at a fixed mass. If we take only the four less mas- 
sive objects we find a quite steep relation between metal- 
licity gradient and galaxy mass, parallel to the locus of the 
galaxies with the steepest gradients (we call it the maxi- 
mum steepness boundary line) similar to the predictions of 
the earlier monolithic collapse models. This finding seem to 
be in qualitative agreement with the observational results 
by Spolaor et al. (2009). As for the points near the max- 
imum steepness boundary, they always refer to the models 
with the highest SF efficiency at that given mass. We note 
another trend, symmetric to maximum steepness boundary 
with respect to the solid line (trend of the entire sample) , in 
the sense that at the highest masses we have also the flat- 
test gradients. This seems to go in the direction of Ogando 
et al. (2005), Spolaor et al. (2009) and Jablonka et al. (2007) 
results. In particular, the scatter is minimum at masses be- 
low ~ 10^° Mq. These galaxies tend to have neither shallow 
metallicity gradients nor very steep ones. 

In order to explain such findings, we first note that 
the formal linear regression to our model predictions gives 
^z/H ~ —O.Oiloga, namely a value much smaller (in ab- 
solute value) than the slope of the mass-metallicity rela- 
tion [Z/H]core ~ 0.57 log a. Therefore the relation between 
gradient slope and galactic mass cannot be explained by 
the Carlberg (1984)'s argument (c.f. Introduction, see also 
Jablonka et al. 2007). In other words, the steepening of 
the gradient with mass is not due to the sole increasing 
metallicity of the galactic core, whereas the outermost re- 
gions of galaxies differing in mass keep the same value 
for [< Z/H >]eff. Indeed it has been shown observation- 
ally that the metallicity of the entire galaxy should obey 
to the mass-metallicity relation (e.g. Graves et al. 2007). 
Such a relation is satisfied by our models, for which we pre- 
dict [Z/H]eff ~ 0.53loga. Hence Az/h "^ l^/H]eff(o') — 
[Z/H]corei(T) ~ 0.53 log a - 0.57 log a = -0.04 log a'^ The 
reason for this increase in the global galaxy metallicity with 
mass is due to the fact that the entire galaxies, not only 
their central cores, should form more efficiently as their 
mass increases in order to comply with the downsizing trend, 
namely the need to have [< a/Fe>] ratios greater than zero 
and positively correlated to the mass. This request renders 
the average gradient slope predicted by the revised mono- 
lithic models flatter than the earlier monolithic collapse 
models a la Larson. However, galaxies with steep gradients 
still exists (e.g. models La and Lb) and lie on the maxi- 
mum steepness boundary. On average, galaxies with mass 
~ IO^^Mq feature gradient slopes quite close to the maxi- 
mum steepness boundary, therefore the scatter is small. At 



^ Note that in our simulations log{Rcore / Ref f) 
ity of the cases. 
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larger masses, the average gradient is nearly one half of the 
maximum steepness boundary value at that mass, hence al- 
lowing for more intermediate possibilities. 

It is interesting to understand what are the major 
causes for such a range although we have unevenly sampled 
the parameter space and despite the not very high number 
of simulated galaxies. We suggest that differences in the ini- 
tial conditions of the protogalactic cloud(s) can reproduce 
the observed scatter. Here we discuss their relative role. 

• Especially at larger masses, the higher the star forma- 
tion parameter esp, the steeper the gradient (see Fig. [S]). 
As an example compare the model Ma2 with Ma3 (or mod- 
els E6 and E2c). Their initial conditions are the same but 
for tsF- The metallicity gradient predicted for the former 
model (which has a higher star formation parameter) is 
much steeper than the latter case. At the first order, the 
metal production scales with the SF rate. Therefore, the 
metallicity increases faster when esp is increased and all the 
other parameters are held fixed. However, the SF depends 
also on the cooling, that increases at higher metallicity, and 
the local gas density. Another important factor in regulating 
the SF is the interplay between stellar feedback and local po- 
tential. Taking togheter all these factors, in the case in which 
the gas is already in place, the net product is that augment- 
ing esp leads to a faster increase (relatively speaking) in the 
metallicity of the central regions, where the potential is very 
deep and the gas is denser than in the outskirts, where an 
higher SF rate also implies that the conditions for the wind 
set in earlier. 

• The role of the chosen profile is slightly less evident. 
On average, the flat profile leads to slightly steeper gradi- 
ents. For instance, compare model Mb4 with Mai that share 
the same initial condition but the profile. Similarly, compare 
model Elb to Ela. This happens because, while in the IS 
models most of the (pristine) gas is already in place, in the 
flat models the majority of the gas supply to build the inner 
regions has go through the outskirts when being accreted. 
Hence, the sinking gas is polluted by metals, leading to a 
faster metal enrichment of the inner regions. However, this 
effect is weaker than that caused by esp- For instance com- 
pare model Ma2 with Mb4. 

• As for the temperature, starting from a higher value 
implies a longer time for cooling the gas and feeding the 
star formation process. In a sense, the effect is similar to 
the difference between the flat case versus the IS case. For 
instance, on the basis of the previous point we would ex- 
pect model Mai to exhibit a (slightly) shallower gradient 
than the one of model Ma3. Instead it is steeper. However, 
a higher initial temperature is not enough to counterbalance 
the effect of a large change in esp (see model Mai versus 
model Ma2). 

• For flat models, the initial gas density seems to be rel- 
atively unimportant (e.g. compare models E2c and E4c) in 
the determining the slope of the metallicity gradient. 

In conclusion, we do not find a parameter that fully governs 
the creation of the gradient, even if esp seems to be quite 
important. Different - but reasonable - combination of the 
input parameters lead to model properties that obey both 
the overall properties observed in elliptical galaxies and ex- 
hibit average metallicity gradient of -0.3 dex per decade in 
radius. Changes in the initial conditions within the same 



broad formation scenario create the scatter in the predicted 
gradients at a single mass. These changes, therefore, should 
not be ascribed to different pictures for the formation of 
the galaxies. They rather mimic cases in which the accre- 
tion from the proto-galactic clouds may be faster (e.g. the 
IS cases) or proceed through cold accretion through fila- 
ments (Dekel & Birnboim, 2006, the flat case). They also 
show the different behaviour of models where the star for- 
mation is favored (higher esp, e.g. for the formation of the 
most massive galaxies) or disfavored (models with high ini- 
tial temperature: the gas is accreted in pre-existing haloes 
and has to cool before forming stars). 

In the middle and right panels in the upper row of Fig. |4] 
we compare our model predictions to metallicity gradients 
inferred from observations (Mehlert et al 2003, Annibali et 
al. 2007, Jablonka et al. 2007). Obviously the above discus- 
sion on the cause of the (scatter in the) metallicity gradient 
applies also to the other mass tracers (cr and the central 
[< O/Fe >]). As explained above, however, here we can 
compare our predictions with the values measured by the 
observers. We can thus show that the predicted range as 
well as the average gradient slope (-0.3 dex per decade in 
radius) are in agreement with observations. We note how 
different observational groups infer slightly different mean 
/S.Z/H (e-g- compare the samples in Fig. |4]). For instance in 
the literature average values either as low as —0.22 ± 0.1 
or as high as -0.34 ± 0.08 (Brough et al., 2007) can be 
founqj. Still consistent with each other, though. This might 
be due to a different combination of line-strength indices 
used to infer the variation in metallicity (see the analysis in 
Sanchez-Blazquez et al. 2006) . Also differences in the SSP li- 
brary used to transform indices into abundances can create 
the offset. Moreover, small number statistics can still bias 
the results as well as the fact that, even in the same sam- 
ple, metallicity gradients are not measured out to the same 
radius. Some authors claim the difference is caused by the 
environment, with field ellipticals featuring shallower gra- 
dients on average with respect to galaxies living in higher 
density regions (Sanchez-Blazquez et al. 2006). Such a sug- 
gestion might explain the offset between the Mehlert et al. 
(2003) Coma cluster ellipticals and the Annibali et al. (2007) 
spheroids. 



^ We refer to Table 4 in Spolaor et al. (2008, and references 
therein) for a useful comparison of the gradients in age, metallic- 
ity and a- enhancement inferred by the above mentioned obser- 
vations. 
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Figure 4. Theoretical predictions versus abundance gradients 
inferred from observations. Upper panels: Predicted [< Z/ H >] 
gradient slope in the stellar component of our model ellipticals 
(hollow circles) and bulges (full dots) as a function of the stellar 
mass and other mass tracers (namely the stellar velocity disper- 
sion and the central [< O/Fe >]). Formal linear regression fits 
to the entire sample (solid line), to the galaxies with steepest 
gradients (dotted line, the maximum steepness boundary) and to 
dwarf ellipticals (dashed line) are shown. Data from Mehlert et 
al. (2003, asterisks, red), Annibali et al. (2007, crosses, blue) and 
Jablonka et al. (2007, stars, green) are shown. Lower panels: As 
above, but for the [< O/Fe >] gradient. 
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5.1.2 Gradients in abundance ratios 

We now move to the analysis of the bottom row of Fig. [l] 
No clear relation with mass is found for the [<0/Fe>] radial 
gradients. Indeed, as expected from Paper I and II, most of 
our models predict a nearly flat [< O/Fe >] gradient, with 
some showing either positive or negative slopes. 

In particular, we suggest the gradient in the [< q/Fo] 
ratio to be related to the interplay between the velocity of 
the radial flows moving from the outer to the inner galactic 
regions, and the intensity and duration of the SF formation 
process at any radius. Clearly a larger or smaller parameter 
of star formation can have a strong influence on this process. 
This result implies that we do not need the merger events 
in order to have a shallow [< a/Fe >] gradient. 

In general, we find that in our models with Ao/f,^ < 
the role of both the gas flowing inward and the star for- 
mation timescale increasing at large radii is non negligible. 
The role of the initial temperature can be important. If the 
galaxy formation process starts from hot gas (i.e. 10*^"^ K), 
we predict I^q/Fe > in the majority of the cases. They 
are thus similar to the quasi-monolithic chemical evolution 
models of Pipino & Matteucci (2004, PM04) with not inter- 
acting shells in which the infall timescale increases at shorter 
radii, whereas the star formation efficiency is constant. On 
the other hand, models starting with cold (i.e. 10*~* K) gas 
seem to prefer a negative Aq/fc- 

The sole SF efficiency seems to affect the predicted ab- 
solute value of the gradient slope; in fact all the models 
most effective in forming stars exhibit the steepest slopes 
at the same time. Basically, an increase in the SF efficiency 
enhances the differences between the inner core and the out- 
skirts set by the other initial conditions. For instance, if the 
gas is already in place, a high efficiency in forming stars 
boosts the outside-in process. In such a case, the star for- 
mation process, which also locks the metals into the stars, 
is fast enough in the central regions to avoid the contam- 
ination of the metals flowing from larger radii. In practice 
we end up in the extreme case in which the gas flows can 
be neglected and Ao/Fe ~ 0.2 as in the standard chemical 
evolution models (Pipino et al. 2006). 



5.1.3 Correlations between gradients in metallicity and 
gradients in abundance ratios 

The flnal part of the theoretical analysis involves the study 
of possible correlations between gradients in metallicity and 
gradients in abundance ratios. As a conflrmation of what 
said in Sec. 5.1.2, galaxies showing the steepest positive 
[<0/Fe>] gradient slopes, have also quite a strong radial 
decrease in the [<Fe/II>] ratio (Fig. [S|. These galaxies are 
also the most massive ones. A correlation in this sense seems 
to be conflrmed by the Annibali et al. (2007) data, whereas 
Mehlert et al.'s (2003) galaxies exhibit values for Aq/fe 
constant with Afe/h- A quantitative confirmation needs a 
sample statistically richer. Perhaps more interestingly, nei- 
ther the observations nor the models cover the region with 
^o/Fe < and Afs/h < —0.4: galaxies with the steepest 
metallicity gradients undergo a strong outside-in formation 
process. In galaxies with Aq/fc < - namely, models that 
likely have a local star formation efflciency decreasing with 
galactocentric radius - the stellar feedback is more effective 
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Figure 5. Predicted [< Z/H >] gradient slope as a function of 
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Figure 6. Predicted [< O/Fe >] gradient slope as a function of 
the [< Fe/H >] one. Symbols as in Fig. [i] 



in contrasting the metal-enhanced flows; therefore, the fl- 
nal Afe/h is smaller (in absolute value), hence closer to 
the expectations from models which do not take into ac- 
count gas flows within the galaxy. At the same time we 
predict a paucity of galaxies in the region Aq/fc > and 
^Fe/H > —0.2. More observations are needed to confirm 
this suggestion. A lack of galaxies can also be noticed in the 
upper left corners on the left panels in Fig. |4l 
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5.2 Gradients in line-strength indices 

Before discussing the implication of our results, it is use- 
ful to re-cast the results of Fig. 4 in terms of their obser- 
vational counterparts in order to enable a ready compari- 
son between our model predictions and future observational 
samples. This is done in Fig. [T] 

We transform the predicted metallicity and abundance 
ratios into line-strength indices by means of Thomas et al. 
(2003) SSPs. In practice, we interpolate the Thomas et 
al.'s theoretical library in order to get a value for the indices 
Mgb and Mg2 for each combination of age, metallicity and a- 
enhancement at any given radius. For simplicity we assume a 
fixed 12 Gyr old population and that the age radial gradient 
are always negligible, since our models always predict age 
differences lower than 0.5-1 Gyr. 

In the upper panels of Fig. [7] we present the results for 
the index Mg2 as a function of the stellar mass and stel- 
lar velocity dispersion. Instead of the central [< O/Fe >], 
here we use the central value for the index, since it corre- 
lates with mass. The solid lines are formal linear regres- 
sion fits to the models. The hatched area represents the 
portion of the plane AM(;2-mass (-cr) covered by the data 
from Ogando et al. (2005). First, we note that the trends 
in the theoretical metallicity gradients with mass is con- 
firmed when translated into observables. In particular we 
confirm the presence of a boundary that corresponds very 
well findings by Ogando et al. (2005), whereas the mean 
predicted slope (rfMg2/log(Rcft,*/Rcorc,«) ~ —0.06 mag per 
decade in radius) is shallower than this limit at any mass. 
We recall that a one-to-one correspondence between index 
gradient and metallicity gradient cannot be done, since the 
transformation depends also on [<0/Fe>] and age. This is 
the main reason why we based our interpretation on Fig. [4] 
rather than on Fig. [T] 

As expected, the region of the planes index-mass tracer 
covered by our model predictions overlaps with those of the 
actual data from the observational samples (Mehlert et al 
2003 - red triangles; Annibali et al. 2007 - blue triangles; 
Jablonka et al. 2007, yellow; see also Sanchez- Blazquez et 
al., 2006, not shown here) that we employed in the previous 
sections. 

No trend appears when we show the predictions regard- 
ing the Mpb index (Fig. [T] lower panels). We show this 
discrepancy as a warning: trends can be strongly index- 
dependent (see also Sanchez-Blazquez et al. 2006). The 
same relation between a theoretical metallicity profile and 
the galactic mass can lead to a different gradient-mass be- 
haviour in the observer plane, depending on the chosen in- 
dex. According to Jablonka et al. (2007), such a dichotomy 
between these two indices should be related to Mg is due 
to the C abundance that plays a role in the index strength. 
Also the role of the age gradients cannot be neglected. 



6 DISCUSSION 

In this paper we showed that difference in the degrees of 
dissipation, in the times at which the galactic wind occurs 
and star formation histories alone can explain the observed 
scatter within a quasi-monolithic assembly. At variance with 
other authors (Kobayashi 2004) we do not need differing 



channels (i.e. "truly monolithic" galaxies, "truly hierarchi- 
cal" galaxy and a mixture of these two) to cover the range of 
observed gradient slopes. In a companion paper (Di Matteo 
et al., 2009) we show, instead, that equal mass dry- mergers 
between ellipticals systematically lower (by a factor of ~ 2) 
the slope of the pre-existing gradient. Therefore we argue 
that if one wants to explain the scatter observed by Ogando 
et al. (2005), Spolaor et al. (2009) and Jablonka et al. (2007) 
in the gradient slopes at high masses with the effects of dry- 
mergers one can accommodate only a few of such episodes, 
otherwise we would observe only galaxies with flat gradi- 
ents. This is true unless there is a channel that continuously 
provides galaxies with the steepest gradients (i.e. -0.5 dex 
per decade in radius) that then can undergo mergers. These 
ellipticals clearly cannot come from mergers, otherwise we 
would need progenitors with slopes even steeper than the 
early monolithic collapse models of Larson and Carlberg 
(i.e. -0.5 to -1 dex per decade in radius, see Introduction). 
Hence the majority of ellipticals have presumably formed 
in a monolithic fashion even if we allow some dry-mergers 
to occur. Therefore we refer to dry mergers as possible (but 
not necessary) episodes in the galaxy life which may change 
the gradient rather than to a well-defined channel for galaxy 
formation which co-exists with the monolithic channel as in 
Kobayashi (2004). Similar constraints on the number of dry 
mergers can be obtained by the [< a/Fe>]-mass relation 
(Pipino & Matteucci, 2008). 

As far as the wet merger^j are concerned, it is argued 
that they may steepen the gradients if star formation takes 
place in the metal rich gas tunneled towards the galactic 
core (Hopkins et al., 2008). Such a mechanism, however, 
creates strong features in the age profiles - at variance with 
observations - that may disappear only after several Gyr if 
the galaxy evolves in isolation since then, at the expenses of 
a flattening of the metallicity gradient. Unfortunately, Hop- 
kins et al's simulation are not done in a cosmological context 
and start from very simplistic assumptions (zero metallic- 
ity disks), therefore it is not clear what happens when the 
simulated galaxy undergoes several mergers as predicted in 
the hierarchical formation scenario. According to Kobayashi 
(2004), the steepening of gradients by the secondary star for- 
mation (i.e., wet merger) seems to occur only rarely. More- 
over, we expect these gradients to be erased by subsequent 
dry mergers. Furthermore, it has been shown (Pipino et al., 
2009a) that the hierarchical formation is still incompatible 
with the observed [a/Fe]-mass relation in ellipticals. On the 
other hand, a clear forecast of our model is non-evolving 
metallicity gradients in time, apart from the effects of the 
passive luminosity evolution on the galaxy spectrum. 

Another prediction of our models is the correlation be- 
tween [<0/Fe>] and [<Fe/H>] gradients, in the sense that 
we expect galaxies with the steepest [<Fe/H>] gradients 
to have a very low [<0/Fe>] abundance ratio in the core; 
therefore such galaxies must exhibit a very steep and posi- 
tive [<0/Fe>] gradient. This does not translate into a clear 
correlation between the [<0/Fe >] gradient and the mass 
because of the large scatter which erase any clear signal. 



^ Wet mergers, as opposed to dry mergers, are tliose where gas 
is involved and SF triggered. 
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Before concluding, we wish to discuss some assump- 
tions, limitations and implications of the present study. 

The initial conditions where chosen in order to repro- 
duce the typical present-day colours, SF rates (as observed in 
high redshift progenitors) and central [<0/Fe>], [<Z/H>] 
values for elliptical galaxies. Preliminary exploration of the 
parameter space led us to restrict our analysis to a range 
0.1-10 in the star formation parameter esF, 10 ~ K in the 
initial temperature, 10% SN efficiency. The majority of the 
resulting models show metallicity profiles linearly decreas- 
ing with log radius, stellar mass to light ratios and radii in 
agreement with observations. 

As a first approximation, our bulge models do not take 
into account the presence of a disk. This is justified by 
the following two reasons: the disk forms on a much longer 
timescale (e.g. Zoccali et al., 2006 Matteucci & Brocato, 
1990, Ballero et al.,2007) and current observational samples 
have been derived by accurately selecting edge-on galaxies 
where the contamination from the discs should be minimal. 
However we stress that our results for the bulges might be 
less constrained and robust that those concerning elliptical 
galaxies and further investigation is needed. 

The number of modelled galax;ies is small, and hence can 
suffer from the same small number statistics that bias the 
observations. We therefore avoided any specific prediction on 
the mean trend of both the metallicity and the [<0/Fe>] 
gradient with mass. The formal linear regressions shown in 
the figures are for the mere purpose of guiding the eye and 
the exact positioning of the maximum steepness boundary 
might depend on the portion of the parameter space ex- 
plored. We stress that the main purpose of the paper is to 
show that even in the monolithic formation collapse a range 
in the predicted gradients consistent with observations must 
be expected and that a typical metallicity decrease of 0.2 dex 
out to 1 Reff can be easily reproduced by recent monolithic 
collapse models. This result is robust, because even if the 
models presented here do not cover the entire parameter 
space, the range in predicted gradients cannot be decreased 
by adding other models. The statistical scatter may change; 
however observational samples likely suffer from the same 
small number statistics problem: this is why we avoided 
any detailed statistical analysis in the present study. The 
presence of elliptical galaxies with fiatter gradients than the 
original models by Carlberg (1984) should not be used as an 
evidence for mergers. 

The metallicity profiles of some of the galaxies with the 
steepest gradients slightly depart from linearity (c.f. Fig. 
3). We would predict a shallower gradient than the one re- 
ported in Table 2 if we limited our analysis to the inner 1/3 
Reff, vice- versa the predicted gradient would be steeper if 
we were to consider only the region 1/3 - IReff- There- 
fore, we caution the reader that the conclusions about the 
steepest gradients in our model galaxies and their relation 
with the monolithic boundary might depend on the chosen 
radius. A detailed comparison between model profiles and 
single well studied galaxies over a large mass range will al- 
low us to study the metallicity gradients in their finer details 
and better constrain the models presented here. 

Moreover, while most of our models obey to the mass- 
size relation for ellipticals (e.g. Shankar et al., 2010), some 
galaxies with similar meiss (e.g. compare models Ma2 and 
La) have quite different radii. The former model has a ra- 



dius consistent with those for normal ellipticals of that mass 
(e.g. Shankar et al. 2010), the latter is more typical of an 
early type Brightest Cluster Galaxy (BCG, e.g. Graham et 
al., 1996). We chose not to make any distinction between 
BCGs and normal ellipticals in our models since gradients 
measured in BCGs have traditionally been included in the 
sample as the ones that we use and because there is no 
difference as far as the chemical properties are concerned 
(e.g. von der Linden et al., 2007, Brough et al., 2007). How- 
ever the reader should keep in mind that a structural dif- 
ference between BCGs and normal ellipticals seems to exist, 
and BCGs seem to harbor steep gradients (e.g. Brough et 
al., 2007). Therefore, in light of the special role of BCGs 
(e.g. Pipino et al., 2009b, von der Linden et al., 2007, and 
references therein), further and dedicated observations and 
modelling are required to ascertain if there is any systematic 
difference in the metallicity gradients with respect to more 
ordinary ellipticals and what is the cause. 

Also, we remind that the majority of the observational 
works use Mg as a proxy for the a elements, as can be eas- 
ily observed in absorption in the optical bands giving rise 
to the well known Mg2 and Algt Lick indices. However, the 
state-of-the-art SSPs libraries (Thomas et al. 2003, Lee & 
Worthey, 2009) are computed as functions of the total a- 
enhancement and of the total metallicity. This is true also 
for the stellar tracks, where the O abundance dominates the 
opacity, and hence the stellar evolution. The latest obser- 
vational results (Mehlert et al. 2003, Annibali et al. 2006 
and Sanchez-Blazquez et al. 2007) that we contrasted to 
our predictions in this study, have been translated into the- 
oretical ones by means of these SSPs; therefore the above 
authors provide us with radial gradients in [a/Fe], instead 
of [Mg/Fe]. This is why in this paper we focus on the theo- 
retical evolution of the a elements by using O that is by far 
the most important. 

Here we briefly recall that both O and Mg come from 
the hydrostatic burnings in massive stars, therefore they are 
produced in lock-step. It has been suggested recently (e.g. 
McWilliam et al., 2007) that this might not be true at solar 
(and above solar) metallicities. While this is an important 
effect in detailed chemical evolution studies it has no impor- 
tance when the luminosity weighted properties of a compos- 
ite stellar population are concerned. This happens because 
luminosity averages weigh more the stellar populations at 
lower metallicities (lower M/L), where the differences be- 
tween O and Mg production are negligible (if any). Finally, 
even if the abundance of O and Mg are offset by some fixed 
quantity (i.e. [0/H] = [Mg/H]+const), the predicted gradient 
would be the same. 



7 CONCLUSIONS 

In this paper we study the formation and evolution of el- 
lipticals and bulges by means of a hydrodynamical model 
(c.f. Paper I and II, respectively) in order to understand the 
origin of the observed scatter in the abundance gradients of 
early type galaxies. Here we summarise our main results. 

• We find Az/H in the range -0.5 - -0.2 dex per decade in 
radius with a mean value of -0.3 dex per decade in radius, in 
agreement with the observations (e.g. Kobayashi & Arimoto, 
1999). 



14 A. Pipino et al. 

• In agreement with Ogando et al. (2005) and Jablonka 
et al. (2007), we find that that the scatter in the gradient 
slopes increases as a function of mass. We reproduce such 
a scatter in the observations by means of variation in the 
initial conditions in galaxy models. 

• Model galaxies which behave as the earlier monolithic 
collapse models by Larson (1974) and Carlberg (1984) de- 
fine a maximum steepness boundary in the metallicity (and 
index) gradient slope-mass plane. These galaxies are prefer- 
entially those with the highest star formation efficiency at 
that given mass. 

• No galaxies with gradients steeper (i.e. more negative) 
than the value given by the our predicted theoretical bound- 
ary are observed (Ogando et al., 2005, and Spolaor et al., 
2009, for ellipticals and Jablonka et al., 2007, for bulges). 

• No correlation between Aq/^,. and other galactic prop- 
erties are found, in agreement observations for ellipticals 
(Mehlert et al., 2003, Annibah et al., 2007) and bulges 
(Jablonka et al., 2007). 

• The abundance gradients in the abundances, 
once transformed into line-strength indices, lead 
to dMg2/log{Rcore,t/ Reff,*) ~ —0.06 mag and 
dLogMgb/log{Rcore,*/Reff,*) ~ —0.1 per decade in 
radius, again in agreement with the typical mean values 
measured for ellipticals and bulges. 

• We note that the behaviour of the gradient slope as a 
function of the galactic mass strongly depends on the par- 
ticular line-strength index used. In fact, the predicted Mg2 
index gradient seems to correlate with mass, whereas the 
Mgb index gradient does not. 

• In Paper I we demonstrated that the differential oc- 
currence of galactic winds (outside-in formation) alone can 
explain the existence of the metallicity gradients discussed 
in this paper. Here we add that this mechanism predicts a 
tight correlation between line-strength index and escape ve- 
locity gradients which has been confirmed by recent data 
(see Scott et al., 2009). 

Larger, homogeneous and statistically meaningful ob- 
servational sample of gradients in elliptical galaxies out to 
one effective radius can confirm such a prediction and vali- 
date the model. 
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Figure 7. Comparison between model predictions and data in 
the observer's plane. Upper panels: Predicted Mg2 gradient slope 
as a function of the stellar mass and mass tracers. The hatched 
area represents the portion of the plane AM(;2-mass {-a) allowed 
by the data from Ogando ct al. (2005). Symbols and data as in 
Fig. |4] Lower panels: As in the upper panel, but for the Mg^ 
index. 
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